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In this paper, numerical predictions of acoustic transmission through a 3D stator ob- 
tained using the NASA BASS code are compared with experimentally measured data. The 
influence of vane count and stagger as well as frequency and mode order on the transmis- 
sion loss is investigated. The data-theory comparisons indicate that BASS can predict all 
the important trends observed in the experimental data. 

I. Introduction 

Computational Aeroacoustics (CAA) is focused on the simulation of acoustic fields produced by unsteady 
flows (see Refs. 1-4 for reviews of the progress made and challenges). To achieve this goal, it is necessary to 
accurately simulate the evolution of flow disturbances, which requires the use of high-accuracy spatial and 
temporal numerical schemes. As these schemes have been developed and validated, the application problems 
of interest have also increased in complexity, now typically involving realistic nonlinear flows about complex 
geometries. 

An application of interest to the aircraft industry is the prediction of acoustic transmission through 
turbomachinery blade rows. In the past, analytical theories for unloaded flat plate configurations were 
developed (e.g., Ref. 5), but these do not account for the effect of realistic geometry, loading, and non- 
uniform mean flows. Numerical methods have also been developed for infinitely thin flat plates, 6 and later 
work has incorporated actuator disks to represent the effect of blade loading. 7 These theories have been 
used to validate two-dimensional CAA predictions. 8,9 

As a step towards the ultimate goal of predicting acoustic transmission through realistic fan stages, CAA 
simulations of acoustic transmission through realistic 3D stator geometries with no mean flow have been 
performed. This work relies on a new experimental database acquired at NASA Glenn Research Center 10 to 
serve as a benchmark for assessing the performance of the NASA BASS CAA code. The results highlight the 
ability of CAA to accurately predict acoustic propagation, transmission, reflection, and scattering from 3D 
stator vanes. The cases shown provide verification of the capability of the BASS code to accurately predict 
the effect of varying stator geometries on a range of incident acoustic modes and identify the challenges that 
remain. 
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II. Numerical Solution of the Governing Equations 


The BASS code is an unsteady CAA code, which uses explicit fourth-order time marching schemes 
combined with high order spatial differencing schemes to accurately solve the unsteady nonlinear Euler 
equations. 11-16 

The BASS code solves the nonlinear Euler or Navier-Stokes equations in chain-rule curvilinear form, 
using a finite- difference time-domain approach. For topological flexibility, the BASS code uses structured 
multiblock grids. For portability and maintainability, BASS is written in standard Fortran 2003, and is paral- 
lelized using the MPI-1 17 standard. BASS has been tested on a range of compilers and MPI implementations, 
and has been validated on a number of benchmark aeroacoustics problems. 

The BASS code has a number of spatial differencing schemes, artificial dissipation methods, and explicit 
time marching schemes implemented. For this work, the BASS code was run using the HALE-RK67 time 
marching scheme. 18 An automatic shock-capturing dissipation model 19 was used, based on the explicit filters 
of Kennedy and Carpenter. 20 The Tam and Webb 7-point optimized fourth-order DRP scheme was used 
for the spatial differencing. 21 The nonreflecting boundary conditions of Giles 23 were used at the driver and 
far-held boundaries, combined with a source term used to impose acoustic modes at the driver boundary. 

III. Description of Experimental Test Cases 

Recently, an experimental stator transmission database has been obtained at the NASA Glenn Research 
Center. 10 The database was acquired on the Advanced Noise Control Fan 24,25 (ANCF) test bed that was 
built in the early 1990s 26,27 to evaluate noise reduction concepts and to provide databases for CAA code 
validation. Figure 1 shows the Aero- Acoustic Propulsion Laboratory where the ANCF is located and a photo 
of the ANCF. Only a subset of the database will be described here. 

A Configurable Fan Artificial Noise System 10 (CFANS) was utilized to generate and control circumfer- 
ential modes (m) and radial modes (n) in the absence of a mean how. The modes were measured by the 
Rotating Rake mode measurement system, 28 which was developed at the NASA Glenn Research Center in 
the 1990s to measure turbofan duct acoustic modes. The system is a continuously rotating radial microphone 
rake that is inserted into the inlet of the duct. It provides a complete map of the acoustic duct modes present 
in a ducted fan. The baseline case was a clean duct, with no stator vanes or fan blades. A parametric set 
of modes was generated at the aft driver set, and measured by the rotating rake at the inlet. Fourteen 
stator vanes at a 45° stagger angle or twenty-eight stator vanes at a 20° stagger angle were installed and 
the rotating rake mode measurements were repeated. It is assumed that the source does not change in the 
presence of the stators, so that the acoustic transmission can be computed by subtracting the mode PWL 
measured with stators installed from the mode PWL measured in the clean configuration. Table 1 provides 
the modes generated and their cut-off ratio at the measurement location inside the duct, for the studies 
documented in this paper. The experimental setup is illustrated in Figure 2. 

In the experimental data, there are two primary sources of error. The first source of error is the error 
in the physical measurement, which is within ±1 dB. The second source of error is modeling error, which 
arises from the postprocessing algorithms used to determine the power level of the acoustic modes from the 
measured data. One source of modeling error arises from the presence of acoustic modes propagating in both 
the positive and negative axial direction in the duct. 32 These modes cause modeling error because the single 
row of microphones used in the rotating rake cannot distinguish between the contribution of each of the 
two families of modes; thus, the postprocessed sound power level for a given azimuthal mode is an unknown 
combination of the modes propagating in each direction. 

IV. Computational Details 

As shown in Figure 3, the computational domain has a far-held boundary outside the ANCF inlet, 
continues into the duct and has another boundary at the acoustic driver location inside the duct. The 
incoming acoustic modes are imposed at this boundary. The analytical development of the imposed acoustic 
modes is given in Appendix I. The Giles 23 nonreflecting boundary condition is used to allow outgoing waves 
to exit the computational domain with minimal reflection. 

For all cases, a reference grid is generated using the GridPro/az3000 grid generator. 29 In this reference 
grid, the stators have zero stagger angle for ease of topology definition and grid generation. The grids have a 
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Figure 1. Experimental Facilities at NASA Glenn Research Center: (a) Aero-Acoustic Propulsion Laboratory; 
(b) Advanced Noise Control Fan. 


Table 1. Experimental Measurements at Inlet 


Case 

Mode 

(m,n) 

Frequency 

Cutoff Ratio 
At Inlet 

Cutoff Ratio 
At Driver 

Clean PWL 
(dB) 

14 Vane PWL 
(dB) 

28 Vane PWL 
(dB) 

1 

(-2,0) 

480 Hz 

1.77 

2.01 

111.1 

111.3 

110.7 

2 

(+2,0) 

480 Hz 

1.77 

2.01 

111.1 

111.3 

110.7 

3 

(-6,0) 

960 Hz 

1.44 

1.45 

106.6 

103.6 

106.2 

4 

(+6,0) 

960 Hz 

1.44 

1.45 

106.2 

105.8 

105.5 

5 

(-4,0) 

480 Hz 

1.02 

1.04 

113.0 

108.1 

103.3 

6 

(+4,0) 

480 Hz 

1.02 

1.04 

113.5 

108.4 

103.7 


minimum of 10 grid points per wavelength of the highest-frequency acoustic modes, which is well within the 
resolution capabilities of the DRP scheme. Once the reference grid has been obtained, the AFRL GridWarp 
volume grid deformation tool of Melville 30 was used to rotate the stators to the desired stagger angle. 

V. Results and Discussion 

The BASS code was run for the full set of frequencies and modal inputs for the baseline, 14 vane, and 28 
vane configurations. For each case, a single azimuthal and radial mode was imposed at the driver location, 
with a reference mode pressure amplitude of 1.4 Pa. 

During the computational run, the flow was converged to a periodic state and solution data was saved 
for the final two cycles of the disturbance mode after convergence. This data was then postprocessed using 
temporal and azimuthal Fast Fourier transforms to obtain the radial distribution of the pressure and axial 
velocity perturbations in the duct. A least-squares method was then used to find the amplitudes of all cut-on 
modes, from which the mode power level is calculated. 

A. Evaluation of the Numerical Solution Accuracy 

In the numerical simulations of acoustic transmission, there are a number of possible sources of error. Two 
major sources of error are due to the imposition of the input acoustic modes and spurious reflections of 
outgoing acoustic modes at the driver location. 

For each simulation, the power level of the imposed acoustic mode was checked at a number of axial 
locations in the vicinity of the driver. Even in the cases for which the modes were near cutoff, the power 
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Figure 2. Schematic of the experimental setup for the ANCF stator transmission tests (exhaust rake was not 
used in these tests). 


level of the imposed mode was accurate to within ±1 dB. For the cases with well cut-on modes, the power 
level of the imposed mode(s) was accurate to within ±0.3 dB. 

The errors in the imposed acoustic modes have two sources. First, the actual mode specification may be 
numerically inconsistent, causing the computed mode to have an incorrect phase and/or amplitude. Secondly, 
if there is an outgoing mode at the driver location, a spurious reflection may occur due to the boundary 
conditions; this reflection will modify the amplitude and phase of the imposed mode. 

In order to determine the relative importance of these two sources of error, the solutions were postpro- 
cessed to determine the levels of the modes propagating in the positive and negative axial directions. The 
imposed modes propagate in the axial direction towards the inlet {—x direction), while the outgoing modes 
propagate in the +x direction. These outgoing modes can be generated by the presence of stator vanes, by 
axial variation in the centerbody geometry, and by the inlet geometry. Investigation showed that the errors 
due to spurious reflections from the driver location boundary condition caused the majority of errors in the 
imposed acoustic modes. 

Each solution was tested to ensure that the imposed modes were dominant; investigation showed that 
the solutions had a signal-to- noise ratio of at least 60 dB. In this case, the ’noise’ in the numerical prediction 
is due to truncation error from finite-precision arithmetic. 

B. Comparison of Numerical Results to Experimental Data 

In the experimental database, the power level of the acoustic modes is given at the inlet plane. The power 
level of the acoustic modes at the driver location is not given. 

In the numerical tests, a mode perturbation pressure amplitude of 1.4 Pa was imposed at the driver 
plane. Since the flow perturbations are very small for this simulation, the solution is effectively linear and 
can thus be scaled to compare with the experimental data. 

To accomplish this, the numerically predicted power level of the input mode propagating in the — x direc- 
tion at the inlet measurement location is scaled to match the experimental data for the ’clean’ configuration. 
This gave a numerical prediction for the power level of the input mode at the driver location, which was 
then used to scale the mode power levels identically for all three configurations. Table 2 gives the numerical 
predictions for the mode power levels at the driver location. 
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Figure 3. Schematic of the computational domain for the ANCF stator transmission tests. 


It is important to note that the experimental mode power level has been used to scale the mode prop- 
agating in the — x direction. This inherently assumes that the power level of the mode propagating in the 
+x direction is small in comparison. The numerical predictions were used to quantify the accuracy of this 
assumption. Table 2 shows that this is a good assumption except in the case of a nearly cut-off mode (cases 
5 and 6). 

Finally, in order to obtain a direct comparison for the acoustic transmission losses between the numerical 
and experimental data, the numerical prediction for the input mode power level at the driver location was 
used to calculate the actual experimental acoustic transmission losses. 


Table 2. Numerical Predictions for Clean Configurations 


Case 

Mode 

(m,n) 

Frequency 

Predicted Driver 
Location PWL (dB) 

Predicted Inlet 
Reflection PWL (dB) 

1 

(-2,0) 

480 Hz 

111.1 

90.7 

2 

(+2,0) 

480 Hz 

111.1 

90.6 

3 

(-6,0) 

960 Hz 

106.6 

79.1 

4 

(+6,0) 

960 Hz 

106.2 

78.7 

5 

(-4,0) 

480 Hz 

115.1 

108.7 

6 

(+4,0) 

480 Hz 

115.6 

109.2 


C. Results for Modes (±2, 0) at 480 Hz 

Figure 4 illustrates the propagation of these acoustic modes through the three geometries. The modes have 
a long wavelength compared to the stator gap and are propagating very nearly in the axial direction. For 
this reason, both stator geometries are ’open’ to the modes, which results in very little acoustic reflection 
or transmission loss. Reflection loss is defined here as the difference in the mode power level of the incident 
mode and the reflected mode. 

Table 3 compares the numerical predictions for transmission loss with the transmission loss computed 
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Figure 4. Predicted perturbation pressure contours at the start of a cycle for (± 2 , 0 ) mode at 480 Hz through 
the three ANCF configurations (a: (— 2 , 0 ) clean, b: (— 2 , 0 ) 14 vanes, c: (— 2 , 0 ) 28 vanes d: (+ 2 , 0 ) clean, e: (+ 2 , 0 ) 
14 vanes, f: (+ 2 , 0 ) 28 vanes) 


from the experimental data. For completeness, the numerically computed reflection loss is also included in 
the table. Recall from Tables 1 and 2 that the modes reflected from the inlet have a power level 20 dB below 
that of the incident modes; thus, the measured data should have very low error due to the postprocessing. 

As the comparison shows, the predicted acoustic transmission is within the range of experimental error 
for these test cases. 

Table 3. Predicted and Measured Transmission and Reflection Losses (dB) for Modes (±2, 0) at 480 Hz 


Case 

Mode 

Transmission Loss 

Reflection Loss 

Transmission Loss 

Reflection Loss 


(m,n) 

14 Vane 

14 Vane 

28 Vane 

28 Vane 



(Pred/Exp) 


(Pred/Exp) 


1 

(-2,0) 

-0.2/+0.2 

-17.6 

0.0/-0.4 

-20.1 

2 

(+2,0) 

-0.2/+0.2 

-17.0 

0.0/-0.4 

-19.8 


D. Results for Modes (±6,0) at 960 Hz 

Figure 5 illustrates the propagation of these acoustic modes through the three configurations. As in the 
previous case, these modes have a long wavelength compared to the stator gap, but are propagating at a 
large angle to the axial direction. For this reason, the stator geometries are open to the modes rotating in 
the positive direction and closed to a varying degree to the modes rotating in the negative direction. Thus, 
the direction of rotation has a large impact on the reflection and transmission losses. 

For the 14 vane stator geometry, Tyler and Sofrin theory 33 predicts that predicts that mode energy 
can scatter from the imposed azimuthal mode 6 to the cut-on counter-rotating azimuthal mode 8. Both the 
experimental and numerical data show that this scattering is most pronounced for the (—6,0) mode, but 
much less energetic for the (±6,0) mode, as would be expected. In Figure 5 (b), the interference patterns 
caused by the presence of the counter-rotating mode can be clearly seen. 
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Figure 5. Predicted perturbation pressure contours at the start of a cycle for (±6, 0) mode at 960 Hz through 
the three ANCF configurations (a: (—6,0) clean, b: (—6,0) 14 vanes, c: (—6,0) 28 vanes d: (+6,0) clean, e: (+6,0) 
14 vanes, f: (+6, 0) 28 vanes) 


The numerical results for mode (+8, 0) may be contaminated by the presence of strong spurious reflections 
from the driver location boundary condition. In the numerical data, the mode power level of the scattered 
(+8, 0) mode, propagating in the — x direction from the stator to the driver location, has a mode power level 
of 102.8 dB. This strong mode itself reflects (spuriously) from the driver location boundary, propagating 
back in the -\-x direction towards the stator with an amplitude of 94.4 dB. While there is an unknown level 
of transmission loss as this mode propagates through the stator, it will have a measurable impact on the 
amplitude and the phase of the scattered mode propagating in the +x direction from the stator to the inlet, 
which has a predicted mode power level of 98.3 dB. An improved boundary treatment at the driver location 
is needed in order to reduce or eliminate these spurious reflections. 

Table 4 compares the numerical predictions for transmission and reflection loss with the experimental 
data for transmission loss, as well as the magnitude of the energy in the scattered mode for the 14 vane 
stator. In all cases, the predicted acoustic transmission loss for the imposed azimuthal mode 6 is equal to 
the measured transmission loss, within the range of experimental error. In the case of the scattered counter- 
rotating azimuthal mode 8, the predicted acoustic transmission shows the correct trends for the magnitude of 


Table 4. Predicted and Measured Transmission and Reflection Losses (dB) for Modes (±6, 0) at 960 Hz 


Case 

Mode 

(m,n) 

Transmission Loss 
14 Vane 
(Pred/Exp) 

Reflection Loss 
14 Vane 

Transmission Loss 
28 Vane 
(Pred/Exp) 

Reflection Loss 
28 Vane 

3 

(-6,0) 

-2.9/-3.0 

-19.4 

-1.2/-0.4 

-26.1 


(+8,0) 

-8.3/-9.5 

-3.8 



4 

(+6,0) 

-0.0/-0.4 

-19.3 

-0.1/-0.7 

-26.7 


(-8,0) 

-26.3/-15.2 

-23.8 
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Figure 6. Predicted perturbation pressure contours at the start of a cycle for (± 4 , 0 ) mode at 480 Hz through 
the three ANCF configurations (a: (— 4 , 0 ) clean, b: (— 4 , 0 ) 14 vanes, c: (— 4 , 0 ) 28 vanes d: (+ 4 , 0 ) clean, e: (+ 4 , 0 ) 
14 vanes, f: (+ 4 , 0 ) 28 vanes) 


the scattered mode, with reasonable agreement with experiment. This case highlights the need for improved 
boundary conditions. 

E. Results for Modes (±4, 0) at 480 Hz 

Figure 6 illustrates the propagation of the acoustic modes through the three geometries. The modes are 
very close to cut-off, and propagate at a very large angle to the x axis. As Figure 6 shows, the modes are 
nearly propagating in the azimuthal direction near the inlet. This test case highlights the errors in both the 
measurement and simulation of nearly cut-off modes. 

As Table 2 shows, these modes reflect strongly from the inlet lip. This is expected behavior from nearly 
cut-off modes. However, the presence of strong modes propagating in both the positive and negative x 
directions results in large modeling errors from the experimental postprocessing. Thus, the predicted driver 
power level may have large errors. However, a newly-developed two-microphone rotating rake will allow the 
power levels of modes propagating in the positive and negative x directions to be determined, reducing or 
eliminating this source of error in future experiments. 32 


Table 5. Predicted and Measured Transmission and Reflection Losses (dB) for Modes (±4, 0) at 480 Hz 


Case 

Mode 

Transmission Loss 

Reflection Loss 

Transmission Loss 

Reflection Loss 


(m,n) 

14 Vane 

14 Vane 

28 Vane 

28 Vane 



(Pred/Exp) 


(Pred/Exp) 


5 

(-4,0) 

- 2.1/-6.2 

-23.3 

-1.7/-11.0 

-2.9 

6 

(+4,0) 

-2.1/-6.4 

-23.3 

-1.7/-11.1 

-2.8 
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This case was also very demanding for the numerical simulations, due to the high mode power levels and 
very slow rate of propagation of the imposed modes. In addition, high levels of spurious reflection from the 
driver location boundary were seen. 

Table 5 compares the numerical predictions for transmission and reflection with the experimental data for 
transmission. In these cases, the predicted acoustic transmission does not compare well with the experimental 
data. However, there are high levels of error in both the experimental and numerical data, making a 
meaningful comparison challenging until both the experimental postprocessing and numerical boundary 
conditions are improved. 


VI. Conclusions 

In this work, the NASA Computational Aeroacoustics code BASS was used to predict the transmission 
losses of acoustic modes propagating through a 3D stator. These predictions compared well with the exper- 
imentally measured data. While it must be emphasized that the cases considered did not include mean flow 
through the duct, the varying solidity and stagger angles of the stators made the test cases challenging. 

The results of this study are verification of the ability of the BASS code to accurately predict the 
propagation, transmission, reflection, and scattering of acoustic modes in realistic geometries. With this 
verification completed, this work will continue moving towards the prediction of acoustic transmission for 
blade rows (both rotors and stators) in the presence of realistic mean flows. 
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Appendix I: Description of Incoming Acoustic Modes 

The analysis presented here sketches that given by Envia , 22 and extends it to obtain the density and 
velocity perturbations associated with the incoming modes, which are required for the CAA analysis. This 
development is purposely left in dimensional form, to ensure that the flow variables can be easily nondimen- 
sionalized for use in the CAA code. 

In this development, it is assumed that the flow perturbations associated with the acoustic modes in- 
troduced at the inflow boundary are low amplitude and is thus governed by the polar coordinate linearized 
Euler equations for a uniform axial mean flow: 


v 1 

Pt + Up x + p ( — + (y r ) r # - (vq ) 6 + u a 


= 0 


u t + Uu x + \p x = 0 

P 

(v r ) t + U (v r ) P ~Pr = 0 


(1) 


(ve) t + U (v e ) x + p 0 = 0 


pr 

y 1 

Pt + Up x + vp ( — + (y r ) r + - (ye) e + u a 


= 0 


These equations can be rearranged to obtain the convective wave equation: 


-77 vt2 _ 2 ( 1 1 


Ptt T 2 Up xt T U Pxx it I p r T p rr T Poo T Pxx ) — 0 


( 2 ) 


The perturbation pressure is assumed to be simple harmonic in time and in the azimuthal direction, and 
the product of the axial and radial functions: 


p(x,r,0,t ) = A p f m (r)g m {x)e l(me wt) 
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Substituting in, the convective wave equation becomes: 


((h 2 - « 2 ) 4^ - 2 i"U d -^ - ^ g m ) f n 


—a 


772 <Urr 


dr 2 


1 1 dfrn rri 

' r dr 


— \ 

r 2 J m J 


9m 


A p e i(<md ~ 


A separation constant k 2 is introduced: 


(4) 


(((,» - «?) - XvUfy - - «V) gm ) U t) 

-® 2 + r + (^ 2 “ 7*~) /™) 


9m 


(5) 


and the equation simplifies to: 


((1 - M 2 ) + ((|) 2 - * 2 ) 9m ) /„ \ 

+ + rlJT + ( k2 ~ 7^) fmjf/m J 

where M = = is the Mach number. 

a 

Defining: 


0 


(6) 


fm (r) = f mn («t) 

gives: 


d fm . 1 dfrn ( 2 ™ \ _p 2 d fmn . ^ df rnn , f 2 ™ 

dr 2 r dr \ r 2 ) 171 K d(nr) 2 r d(ftr) \ r 2 

Setting this equal to zero gives Bessel’s equation: 


(fix) 


2 d 2 f n 


+ (Kr) Jfey + ( (Kr)2 _ m2 ) /mn = 0 


d {nrY 

which has a solution of the form: 


d (fix) 


(7) 

(8) 

(9) 


fmn (kc) — Cl J n (^") C 2 Y n (nr ) (10) 

where n denotes the radial mode given by the solution. 

At the inner and outer walls, the radial derivative of f mn is zero (which specifies no flow through the 
wall) : 


ciJ' n C nr hub ) + c 2 Y; ( nr hub ) = 0 (11) 

ciJ'n (. Kr tip ) + c 2 Yf (nr tip) = 0 

which gives the eigenvalue equation for radial mode n that the separation constant must satisfy as: 


Jn (A mn ru P m K mn^hub ) Jji (A mn^hub ) Y n (fcrnn ^tip) — 0 


(12) 


Here, fi; mn is the separation constant that satisfies the eigenvalue equation for azimuthal mode m and 
radial mode n. 
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Combining the constants and applying the boundary conditions, the function f mn is given by Bessel 
functions of the form: 


fmn (fcmn'r') C mn | J rnn (p^mnY) 

The function g mn is given by: 


Jrnn (%n r kfe) A v / \ 

Y^^mnTHu,)) ( m " ' 


(13) 


9mn (x) = e ikmnX 

where the axial wavenumber k mn satisfies the equation: 


(14) 


(1 - M 2 ) k 2 mn + 2 k mn M- + [n 2 mn - (-) 
which has the upstream and downstream running solutions: 


= 0 


(15) 


Moo , / uo [ 


1 - M 2 J \ a V a 
At this point, the perturbation variables can be written as: 


~ (i_ M 2 )K : 


2 

mn 


(16) 


n — A f (k r i piimd + kmnX-ut) 

jj — sipjmn \ hj mn 1 J c 


P = 


VO 


Ap | _2 ) fmn {^mnV) V 


i{mO-\-k rnri x— UJt) 


u = A 


kr 


■P 


= A„ 


P (tv k mn U) 
i'K'mn 


If ( x r \ A( mQ+kmnX-ut ) 

Yj\ I J mn y^rnn 1 J c 


P (tv k mn U) 


f' ( K r \ JimO+kmnX-UJt) 

TT\ j mn yr^mn 1 J e 


m 


- 


P (tv k mn U) v 


If ( K r \ p i{mO+kmnX-UJt ) 

Yj\ I J mn y^rnn 1 J c 


(17) 


These perturbation variables are imposed as an incoming wave at the boundary of the computational 
domain. 

Note that the radial distribution of the perturbation pressure is identical for both upstream and down- 
stream running acoustic modes. 
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